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摘 要 : 空间 VLBI 是 未 来 重要 的 发 展 方向 ， 但 其 观测 系统 和 轨道 等 不 稳定 因素 会 导致 
VLBI 相 关 处 理 机 所 需 的 预报 时 延 模型 与 实际 时 延 相差 较 大 ， 在 这 种 情况 下 很 难 获得 干涉 条 
纹 。 本 文 提出 了 基于 多 重 网 格 的 射电 源 条 纹 搜 索 算 法 ， 可 以 在 不 使 用 时 延 模型 的 情况 下 ,使 
相关 处 理 机 在 较 大 范围 内 从 原始 数据 中 自动 完成 搜索 干涉 条 纹 。 利 用 滑动 窗口 技术 提高 搜索 
范围 ， 多 重 网 格 技术 提高 搜索 效率 。 实 测 数据 表明 ,该 算法 在 能 够 不 依赖 先 验 预报 时 延 模型 
的 情况 下 成 功 获得 VLBI 射 电源 干涉 条 纹 。 
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0 引 

甚 长 基线 干涉 测量 (VLBI) 是 20 世 纪 60 年 代 后 期 发 展 起 来 的 一 种 射电 干涉 技术 , 具有 高 分 
辨 率 、 高 精度 的 特点 ， 在 天 文 观测 、 大 地 测量 和 深 空 探测 等 领域 得 到 了 广泛 应 用 ”。 其 基本 
工作 方式 是 位 于 基线 两 端的 天 线 接收 射电 源 辐 射出 的 电磁 波 , VLBI 相 关 处 理 机 利用 干涉 原理 
从 实测 数据 中 得 到 时 延 以 及 时 延 变 化 率 , 进一步 获取 射电 源 和 有 关 基 线 的 位 置 等 信息 。 因 此 
快速 、 准 确 地 获得 射电 源 的 干涉 条 纹 ， 对 VLBI 的 数据 处 理 至 关 重 要 。 

相关 处 理 机 是 YLBI 技 术 的 核心 数据 处 理 设备 , 预报 时 延 模型 与 实际 时 延 的 误差 直接 影响 
到 干涉 条 纹 的 处 理 效果 , 因此 处 理 机 对 所 需 预 报时 延 模型 有 很 高 的 精度 要 求 。 但 在 实际 数据 
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处 理 中 ， 地 球 模型 误差 、 大 气 模 型 误差 、 仪 器 误差 以 及 站 坐标 误差 等 都 可 能 导致 时 延 模型 不 


准确 甚至 无 法 引导 处 理 机 正常 工作 。 在 常规 地 面 YLBI 中 ， 由 于 站 坐标 较 稳 定 ， 可 以 通过 预先 
计算 获取 预报 时 延 模 型 与 实际 时 延 之 间 的 误差 , 然后 采用 人 工 反 复 尝试 后 , 获得 准确 的 干涉 
条 纹 , 在 一 次 观察 中 , 通常 不 需要 再 次 进行 搜索 ,但 对 于 空间 VLBI, 观测 系统 和 轨道 不 稳定 ， 
误差 是 随时 变化 的 , 情况 与 地 面 YLBI 有 较 大 区 别 , 因此 无 法 通过 预先 计算 来 对 时 延 模 型 进行 
一 次 性 修正 。 目 前 针对 探测 器 信号 的 无 时 延 模型 条 纹 搜索 研究 较为 成 熟 。 上 海天 文 台 在 这 一 
领域 取得 了 大 量 研究 成 果 , 解决 的 主要 思路 是 利用 探测 器 测控 信号 自身 的 信号 特征 设计 条 纹 
搜索 算法 ”。 但 由 于 射电 源 信号 是 高 斯 白 噪声 信号 ， 也 就 决定 了 无 法 直接 利用 信号 特征 进行 
条 纹 搜 索 。 目 前 在 实际 数据 处 理 中 , 无 时 延 模 型 的 射电 源 条 纹 搜索 往往 采用 半 人 工 的 方式 进 
行 ， 主 要 依赖 于 操作 人 员 的 经 验 ， 这 种 方式 费时 且 费 力 。 对 于 未 来 空间 VLBI 的 应 用 场景 ， 这 
种 常规 的 半 人 工 方式 难以 满足 需求 。 因 此 , 解决 不 依赖 时 延 模型 对 射电 源 条 纹 搜 索 这 一 难题 ， 
对 推动 空间 VLBI 的 应 用 具有 重要 意义 。 

针对 无 时 延 模型 的 射电 源 条 纹 搜 索 ， 本 文 提出 了 基于 多 重 网 格 的 射电 源 条 纹 搜 索 算 法 ， 
能 够 在 较 大 的 时 延 和 时 延 率 平面 内 进行 搜索 , 利用 消 窗 技 术 以 及 多 重 网 格 技术 快速 自动 定位 
到 精确 的 时 延 、 时 延 率 。 


1 VLBI 相关 处 理 原 理 


VLBI 观 测 台 站 接收 到 射电 源 信号 ， 经 过 变频 、 滤 波 和 二 次 变频 后 ， 被 划分 成 多 个 频率 通 
， 经 过 采样 量化 编码 后 记录 在 磁盘 上 或 通过 网 络 传输 到 YLBI 中 心 “。 设 台 站 1 接收 到 的 信 
为 x1， 台 站 2 接收 的 信号 为 x»， 两 全 站 信号 之 间 的 关系 满足 : 


Tr 


= 


道 
5 


X = x,(t — T) 


因为 地 球 自 转 和 射电 源 相 对 运动 的 影响 ，r 实 际 上 是 一 个 随 独 时 间 变 化 的 函数 ， 在 很 短 
的 时 间 At 内 Tt 可 近似 成 一 次 函数 : 


T(t) = T, + t,t 
SEH T, EAS BEN THIER oD, t Ate TH) AR LA ABA 
经 过 下 变频 后 两 台 站 的 输出 信号 分 别 为 : 


s1(t) =x, Ge Pot 
So(t) = xX2(t)eJ2™ot = x(t — Te J2* fot 


其 中 名表 示 天 空 频率 。 
将 原始 数据 解码 后 , 固定 台 站 1 信号 , 对 台 站 2 信号 做 时 延 补偿 , 设 ts 是 搜索 的 补偿 时 延 ， 
它 与 实际 时 延 tj 之 间 为 差 为 Ats， 表 示 为 : 


Ts = Tr — AT, 


由 于 是 以 数字 信和 号 的 形式 存储 数据 ， 采 样 间隔 是 AT， 则 rs 由 整数 部 分 npAT7 和 人 小数 部 分 
Tsf 两 部 分 组 成 : 


Ts = NAT + Tor 
时 延 的 补偿 分 为 两 部 分 , 一 部 分 是 在 时 域 完成 的 整数 比特 补偿 ， 另 一 部 分 是 在 频 域 完成 
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的 小 数 比 特 补偿 。 对 台 站 2 信号 进行 整数 比特 补偿 后 为 ”: 


sz(t — nAT) = x(t — nAT)e ~i? -nAT) 

= x2 (t — Ty —Tsf + ATts)e Jarolt — Tr + Tsf + ATs) (7) 
通过 条 纹 旋转 ， 消 除 时 延 中 时 变 项 ， 引 入 e-jzrptst 因子 ARRINENER 与 实际 时 
延 率 之 间 的 差 At 满足， 


t, =T, —At, 


` 


q: 


进行 储 里 叶 变 换 后 得 到 其 频谱 $5,(f)， 完 成 小 数 比特 补偿 后 
S27s(f) = X(f + fo)e 12M oats @— jam foatsAt 


对 台 站 1 做 傅 里 叶 变 换 ， 得 到 两 台 站 之 间 的 互 功率 谱 为 : 
Piz =< $1(f),52(f) > (10) 
= |X (fre ee ZTfoAtsAto-j(p1+92) 
其 中 pi 和 2 分 别 是 两 台 站 信和 号 的 原始 相位 。 
通过 上 述 分 析 可 以 看 出 ， 两 台 站 互相 关 功 率 谱 和 补偿 时 延 rs 与 实际 时 延 的 差 值 
4ts， 补 偿 时 延 率 和 实际 时 延 率 的 差 值 4 有 关 , 因此 经 过 不 同时 延 和 时 延 率 补偿 后 的 两 台 
站 信号 互相 关 功率 谱 也 就 不 同 。 当 4rs 和 4 如 为 0 时 ，si() 和 sz() 两 台 站 信号 的 互相 关 功 率 
谱 有 最 大 峰值 ， 此 时 经 过 补偿 后 的 两 台 站 信号 为 同一 波 前 信号 。 
2 多 重 网 格 的 条 纹 搜索 
基线 两 端的 天 线 接收 由 射电 源 辐 射出 的 电磁 波 , 因为 地 球 和 射电 源 的 相对 运动 使 得 电磁 
波 的 波 前 到 达 两 个 天 线 的 时 间 差 值 不 断 改 变 , 干 涉 条 纹 只 存在 于 同一 波 前 信号 的 互相 关 功 率 
谱 中 “。 因 此 ， 只 有 将 接收 到 的 两 台 站 信号 进行 准确 的 时 延 、 时 延 率 补偿 才 可 以 得 到 射电 源 
的 干涉 条 纹 。 本 文 提出 一 种 基于 多 重 网 格 技术 的 算法 , 可 在 无 预报 时 延 模 型 情况 下 对 射电 源 
条 纹 进 行 搜索 。 
2.1 多 重 网 格 技术 基本 原理 
多 重 网 格 的 基本 思想 就 是 建立 多 层次 不 同 尺寸 的 网 格 , 通过 这 些 网 格 在 不 同 的 范围 内 搜 
索 我 们 需要 的 目标 值 ”。 如 图 1 所 示 ， 搜 索 空间 被 划分 为 L 层 ， 每 层 的 网 格 都 由 边 长 相同 的 单 
元 格 构成 ， 单 元 格 的 边 长 称 为 搜索 步 长 ， 每 个 单元 格 的 结 点 表示 不 同 的 搜索 值 。 在 划分 过 程 
中 ， 随 着 网 层 号 的 降低 ， 网 格 的 尺寸 减 小 ， 搜 索 精 度 提 高 。 


~ 


All 多 重 网 格 示意 图 
Fig. 1 Schematic diagram of multiple grids 


通过 对 VLBI 相 关 原 理 的 分 析 , 可 知 每 一 对 补偿 时 延 和 补偿 时 延 率 都 会 与 其 互相 关 功 率 谱 
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相对 应 。 条 纹 搜 索 需 要 在 时 延 - 时 延 率 平面 中 将 存在 条 纹 的 互相 关 功 率 谱 所 对 应 的 时 延 和 时 
延 率 找 出 。 将 时 延 -时 延 率 搜索 平面 划分 网 格 ， 网 格 尺寸 的 大 小 也 就 决定 了 时 延 和 时 延 率 的 
搜索 精度 ， 最 底层 的 网 格 尺寸 由 YLBI 相 关 处 理 机 的 精度 需求 决定 。 

VLBI 处 理 机 对 时 延 和 时 延 率 的 补偿 精度 有 很 高 要 求 , 否则 无 法 获得 条 纹 。 时 延 补偿 的 精 


度 必须 满足 ”: 
|Ar| < -Z> (11) 
其 中 ，BW 为 频率 通道 带宽 ,N 为 谱 通道 个 数 〈FX 型 处 理 机 ) 或 时 延 通 道 个 数 〈XF 型 处 理 机 )。 
时 延 率 的 误差 Ai 需要 满足 : 
IAi| < zis (12) 


其 中 ，7 为 积分 时 间 ，f 为 观测 频率 。 

若 直 接 选 取 满 足 相关 处 理 机 精度 要 求 的 单一 网 格 , 相应 的 计算 量 极 大 , 条 纹 搜 索 效率 很 
低 。 为 提高 搜索 效率 , 采用 粗 网 格 定 范围 , 细 网 格 定 终 值 的 策略 。 首先 选 取 较 大 的 搜索 步 长 ， 
在 粗 网 格 的 搜索 平面 内 确定 正确 的 时 延 和 时 延 率 补偿 所 在 范围 , 如 图 2(a), 通过 计算 确定 正 
的 值 落 在 图 中 红色 区 域 位 置 ,用 于 第 二 步 精 搜索 ,这 一 过 程 称 为 粗 搜索 ,缩小 搜索 范围 后 ， 
同时 缩小 网 格 以 满足 相关 处 理 机 精度 要 求 ， 如 图 2(b) 所 示 ， 进 而 可 在 细 网 格 内 搜索 到 更 高 
| 精度 的 时 延 和 时 延 率 补偿 ， 完 成 精 搜索 。 
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(a) 粗 网 格 示 意图 (b) 细 网 格 示意 图 
图 2 网 格 操作 示意 图 
Fig.2 Schematic diagram of grid operation 

射电 源 信 号 的 信 噪 比较 低 ， 需 要 增加 积分 时 间 来 解决 这 一 问题 , 但 计算 机 对 相关 计算 长 

度 和 FFT 计算 的 长 度 都 有 限制 ”"”， 通 过 对 信号 的 分 段 处 理 可 平衡 两 者 之 间 的 关系 。 在 两 路 信 

号 中 ， 需 要 固定 一 路 信号 作为 参照 信号 , 对 另 一 路 信号 进行 分 段 处 理 。 第 一 路 的 信和 号 与 第 二 

路 的 信号 每 段 进行 相关 处 理 即 实现 滑 窗 相 关 。 信号 滑 窗 操作 方法 如 图 3 所 示 , 将 台 站 信号 在 

搜索 范围 内 以 网 格 要 求 的 搜索 步 长 滑动 , 并 在 相关 计算 后 记录 互相 关 功 率 谱 , 提取 相位 后 计 
算 方 差 。 


数据 窗 
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延 。 如 图 4(a) 表示 时 延 和 时 延 率 补偿 与 实际 时 延 和 时 延 率 相近 的 相位 结果 , 若 不 能 补偿 正确 


图 3 滑 窗 操 作 示 意图 
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Fig.3 Schematic diagram of sliding window operation 
2. 2 评估 时 延 和 时 延 率 补偿 结果 的 依据 和 方法 
由 公式 10 可 知 ， 若 补偿 准确 ， 互 相关 功率 谱 的 相位 与 频率 呈 线 性 关系 ， 和 斜率 表示 残余 时 


的 时 延 和 时 延 率 ， 相 频 图 将 会 如 图 4(b) 所 示 。 
从 图 4 可 以 看 出 ， 互 相关 功率 谱 的 相位 弥散 程度 可 以 衡量 时 延 、 时 延 率 补偿 的 准确 性 。 
补偿 时 延 和 时 延 率 与 实际 时 延 和 时 延 率 误差 小 ， 相 位 弥散 程度 小 ; 补偿 误差 大 ， 相 位 弥散 程 


度 大 。 


要 度量 相位 的 弥散 程度 需要 先 对 相 频 图 进行 线性 拟 合 , 定 下 基准 点 的 位 置 。 在 线性 拟 合 


时 ， 采 用 最 小 二 乘法 实现 拟 合 ， 如 图 4(a). A 


日 于 最 小 二 乘 的 曲线 拟 合 方法 是 通过 最 小 误差 


的 平方 找到 拟 合 结果 , 虽然 补偿 不 准确 的 相 频 谱 不 具备 线性 关系 , 但 仍然 可 以 拟 合 出 一 条 与 
真实 数据 距离 最 近 的 直线 ， 如 图 4(b) 。 


时 延 时 延 素 补偿 准确 


phase(rad) 


0.5 1 


(a) 补偿 正确 的 相 频 图 拟 合 结果 


在 相 频 图 中 ， 
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5 (a) 所 示 。 
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Fig.4 Fit 
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图 4 在 相 频 


ting result on 


时 延 时 延 率 补偿 不 准确 


phase(rad) 
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(b) 补偿 不 正确 的 相 频 图 拟 合 结果 
图 上 的 拟 合 结果 


the phase-frequency diagram 


Ul: 相位 卷 绕 。 由 公式 10 ATS, ARR EHO Fe 
得 到 ， 可 得 相位 是 相 加 的 关系 ， 表 示 为 : 
Of) =p) + O2(f) 


go (FEET EZ Ta), (EP AA EE GENE Cr, n) ZA H 
于 计算 机 在 处 理 相位 的 时 候 只 能 在 Cr, n) 中 表示 ， 此 时 相位 图 被 截断 ， 发 生 跳 变 ”。 如 图 
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(a) 相位 卷 绕 (b) 相位 解 卷 绕 
图 5 卷 绕 与 解 卷 绕 
Fig.5 Phase wrapping and phase unwapping_ 

从 图 中 可 以 看 出 条 纹 是 存在 的 ， 但 相位 发 生 了 截断 ， 发 生 条 纹 平行 出 现 的 现象 ， 导 致 无 
法 对 相位 数据 进行 准确 的 线性 拟 合 ， 过 大 的 拟 合 方差 也 会 影响 后 续 的 判断 。 
通过 上 面 的 分 析 知 ， 因 为 两 台 站 信号 互相 关 功 率 谱 的 相位 超出 (Cr, m) 范围 导致 相位 卷 
绕 现 象 。 为 了 可 以 准确 表示 出 互相 关 功 率 谱 的 相位 ， 需 要 对 相位 进行 解 卷 绕 处 理 。- 解 卷 绕 
方法 为 : 规定 跳 变 阔 值 ， 当 相位 跳 变 超过 闪 值 时 ， 将 相位 进行 搬移 ， 使 截断 的 相位 连接 在 一 
起 ”。 解 卷 绕 后 的 相位 图 (b) 所 示 。 


2.3 多 重 网 格 条 纹 搜索 流程 


图 6 是 多 重 网 格 条 纹 搜索 的 具体 流程 。 首 先 需 要 计算 相关 处 理 机 的 精度 要 求 ， 确 定 粗 搜 
索 的 时 延 和 时 延 率 补偿 的 搜索 步 长 ， 读 取 台 站 记录 的 信号 ， 解 码 后 在 绝对 时 间 上 对 齐 ， 固 定 
台 站 1 信号 , 对 台 站 2 信号 的 进行 粗 搜 索 的 滑 窗 处 理 , 完成 时 延 和 时 延 率 的 补偿 后 计算 互相 关 
功率 谱 ， 搜 索 到 互相 关 功 率 谱 弥散 程度 最 小 的 位 置 。 缩 小 搜索 范围 ， 以 满足 精度 要 求 的 搜索 
步 长 进行 相关 处 理 ， 在 信号 相关 后 得 互相 关 功 率 谱 ， 提 取 相 位 。 计 算 每 组 时 延 、 时 延 率 对 应 
相位 的 方差 并 记录 ， 最 后 将 记录 的 方差 结果 绘制 输出 ， 选 取 峰 值 对 应 的 时 延 、 时 延 率 ， 最 终 
得 到 精确 的 时 延 、 时 延 率 并 输出 。 
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mo NRE. IN ER AME 
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a 
图 6 多 重 网 格 条 纹 搜索 流程 图 


Fig.6 Multi-grid fringe search flow chart 
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3 实测 信号 验证 

为 了 验证 多 重 网 格 条 纹 搜索 的 正确 性 ， 我 们 以 嫦娥 四 号 探 月 工程 中 一 段 实测 数据 进行 
验证 ， 其 观测 代码 为 s88b06a， 观 测 源 为 3C273B。 任务 采 用 的 相关 处 理 机 为 FX 型 处 理 机 。 本 次 
实验 选择 天 马 站 和 乌鲁木齐 站 两 台 站 的 信号 。 信和 号 的 采样 率 为 4MHz， 采 用 2bit 量 化 , 划分 为 
16 个 通道 ， 通 道 带宽 为 2MHz， 积 分 时 间 2s。 本 次 实验 选取 第 一 通道 的 信号 ， 中 心 频 率 万 为 
2233. 5MHz。 由 公式 11 和 公式 12 计 算 可 知 ， 时 延 、 时 延 率 补偿 的 精度 需要 满足 : 


16 
|At| < FX IMHZ =4x10-es 


1 
It < 5s x 2333MHz 


算法 本 身 对 时 延 和 时 延 率 的 搜索 范围 在 理论 上 没有 硬性 约束 。 但 为 方便 起 见 , 在 实验 操 
作 时 ， 以 先 验 模 型 为 初 值 确定 了 较 大 的 搜索 范围 。 在 实际 工程 应 用 中 ,， 若 有 更 大 的 搜索 范围 
需求 而 带 来 更 大 的 计算 量 ， 可 通过 MPI、CUDA 等 并 行 计算 技术 实现 解决 。 


首先 选择 粗 网 格 进行 搜索 :时 延 搜索 步 长 为 1 x 10-5s, 时 延 率 搜索 步 长 为 1 x 10-?s/s。 
在 搜索 范围 内 ， 分 别 计算 不 同 的 时 延 和 时 延 率 所 对 应 的 互相 关 功 率 谱 的 相位 方差 c2, 为 了 使 
结果 更 加 直观 ， 取 a? 的 倒数 绘制 成 图 7 (a) 。 得 峰值 位 置 为 补偿 时 延 : 2.88 x 10-3s ， 补 偿 
时 延 率 :7.15 x 10-7s/s。 因 此 时 延 补偿 的 范围 是 2.87 x 10-3s~ 2.89 x 1073s, IN MER AME 
范围 为 : 7.14 x 1077s/s~7.16 x 10-7s/s。 后 面 的 条 纹 搜索 在 此 范围 内 进行 更 为 精确 的 搜索 ， 
得 到 时 延 、 时 延 率 补偿 的 精确 值 。 


~ 1x1071s/s 


时 延 尝 (s/s) mE 24 


时 延 (s) 


(a) 粗 网 格 搜 索 图 (b) 细 网 格 搜索 图 
图 7 多 重 网 格 条 纹 搜索 


Fig.7 Multi-grid fringe search 


经 过 上 述 的 操作 ， 确 定 精 搜索 时 延 精 搜索 范围 为 : 2.87 x 1073s~ 2.89 x 1073s, 时 延 率 
精 搜索 范围 为 7.14 x 10-7s/s~7.16 x 10-7s/s。 精 搜索 时 延 步 长 为 1 x 1076s, -时 延 率 搜 
索 步 长 为 1 x 10-10s/s。 分 别 计 算 不 同时 延 、 时 延 率 对 应 的 Q 值 ， 绘 制 成 图 7(b) ， 求 得 峰值 
对 应 的 时 延 值 是 2.881 x 10-3s, 时 延 率 是 7.151 x 10-7s/s。 将 搜索 到 的 时 延 、 时 延 率 带 入 相 
关 处 理 机 中 进行 验证 ， 得 到 射电 源 的 条 纹 ， 如 图 8 所 示 。 


0.5 1 1.5 2 25 3 3.5 4 4.5 5 
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图 8 实测 数据 条 纹 
Fig.8 Fringe phases of actual data 


本 文 实验 所 使 用 的 的 CPU 型 号 是 Intel(R)Xeon(R)CPU X5650@2.67GHz。 多 重 网 格 条 纹 搜 
索 在 本 次 实测 数据 应 用 中 如 窗口 大 小 ， 程 序 运行 时 间 等 相关 参数 如 表 1 所 示 。 
Rl 多 重 网 格 条 纹 搜索 的 相关 参数 
Tab. 1 Parameters of Model-free Fringe Searches 
参数 名 称 数据 
时 延 搜索 窗 大 小 (s) 


1x103 
时 延 搜 索 精度 (s) mE 
时 延 率 搜索 窗 大 小 (s/s) EE 
时 延 率 搜索 精度 (s/s) mie 
搜索 时 间 (s) 


4 结论 


本 文 针 对 射电 源 VLBI 自 动 条 纹 搜 索 这 一 难题 , 基于 多 重 网 格 的 方法 , 提出 了 一 种 射电 源 
信号 的 条 纹 搜索 算法 , 并 通过 嫦娥 四 号 实测 信号 进行 验证 。 该 方法 解决 了 射电 源 条 纹 搜索 不 
依赖 于 时 延 模 型 这 一 难题 .基于 该 算法 的 应 用 背景 ,经 过 后 续 深 入 的 研究 , 有望 用 于 空间 VLBI 
数据 处 理 中 。 
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Abstract: Space VLBI is an important development direction in the future, but its 
observation system and orbital instability will cause the forecast delay model to differ greatly from 
the actual delay. In such case, it is difficult to obtain interference fringes. This paper proposes a 
radio grid fringe search algorithm based on multiple grids, which can search interference fringes 
from the original data over a wide range without using a delay model. Use sliding window 
technology to increase search range, and multi-grid technology to improve search efficiency. The 
observation data experiment results —show that the algorithm successfully obtains the interference 


fringes of VLBI radio source without forecasting delay model. 


Key words: VLBI; Radio source; Fringe search; Multi-grid; Correlator 
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